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We study the effects of next-nearest neighbor (NNN) interaction on the dynamic phase transition 
(DPT) and hysteresis loop area law in the two-dimensional ferromagnetic kinetic Ising model. We 
find that inclusion of the NNN interaction causes the DPT boundary of the NN kinetic Ising model 
to shift to larger values of magnetic field and temperature. The NNN kinetic Ising model can 
therefore exhibit an interaction induced DPT. Also in the low frequency limit (f^O) the hysteresis 
loop area law, A(ho,f), changes from h^^-^^f'^^ (NN) to hO i^^o oifO i^^o-oi (NNN) where ho is the 
external field amplitude and f is the frequency. DPT and hysteresis in the kinetic Ising model arises 
as a competition between the system's metastable lifetime and the time period of the external field. 
Including the NNN interaction changes the system's metastable lifetime. This causes the DPT and 
the hysteretic properties to change. We conclude that the systems metastable lifetime is sensitive not 
only to the lattice size, external field amplitude, and temperature but also to additional interactions 
present in the system. Furthermore by tracking the probability density P(Q) of the dynamic order 
parameter Q we conclude that the DPT is of second order. 

PACS numbers: 64.60.Ht, 64.60.Qb,75.10.Hk,75.40.Gb,05.40.-a 



I. INTRODUCTION 

The two-dimensional (2D) ferromagnetic nearest- 
neighbor kinetic Ising (NNKI) model has been used ex- 
tensively to study non-equilibrium (NEQ) properties. 
The NNKI model is composed of the Ising Hamiltonian 
and a term representing an external oscillatory mag- 
netic field jl El El H m E]. The presence of an os- 
cillatory field introduces an explicit time dependence in 
the Hamiltonian. This causes the system to exhibit a 
hysteretic response which is a characteristic feature of 
NEQ phenomena. Below the critical temperature (Tc) 
the system exhibits a dynamic phase transition (DPT) 
[i|6l[7l[8l|9l[lQl[ni[Il[ia|14] based upon the strength 
of the field amplitude (ho), frequency (f), and tempera- 
ture (T) of the system. 

Hysteresis in this model occurs due to the competition 
between the system's metastable lifetime and the time 
period of the external field [6 . The hysteresis loop area 
law for the model in 2D, 3D, and 4D has been investi- 
gated using various methods such as Monte Carlo simula- 
tions H El EmSl [16], mean-field theory [H l US [H [IT] 
and 0(N) type models [3]. Within Monte-Carlo sim- 
ulations it has been found that the loop area, A(ho,f) 
oc h^-^f^-^^ in the low frequency limit f^O at a constant 
temperature [5j. The scaling relations for the loop area 
have also been measured in several ultrathin and thin 
ferromagnetic film systems, such as Fe/Au(001) [18], 
FesoNigo [E], Co/Cu(001) [20ii21j, Fe/GaAs(001) [22], 
and Fe/W(110) [l3]. 

Monte Carlo simulations and mean field studies of hys- 
teresis in this model show that the loop area undergoes 
a transition from a symmetric shape to an asymmetric 
shape or vice- versa [HI [16] . This can happen when 
the driving frequency of the field increases (decreases) 



beyond a threshold value or when the external field am- 
plitude decreases (increases) for a particular temperature 
below the critical temperature (Tc) of the model. The 
breakdown in the shape of the hysteresis loop has been 
attributed to a DPT with a spontaneously broken sym- 
metric phase m [5]. A measure of the DPT is given 
by the dynamic order parameter, Q, which is the period 
averaged magnetization. The DPT occurs between an 
ordered dynamic phase with {\Q\) ^ and a disordered 
dynamic phase with [1]. The DPT in the 

NNKI model is of second order [6l [121 US] • The order 
of the phase transition was confirmed by tracking the 
probability density, P(Q), of Q. 

The hysteretic response mechanism and the various 
phases which arise in the NNKI model is dependent upon 
the combination of ho, f, and T values. These parame- 
ters dictate whether the system will be in one of the four 
possible regions: strong field (SF), multi-droplet (MD), 
single-droplet (SD), or co-existence (CE). The SF and 
the MD regions exist for strong fields/large system sizes 
and the SD exists for weak fields/small system sizes [6l. 
Droplet theory as a mechanism for the NEQ response is 
valid in the MD and SD regions but breaks down in the 
SF region [12 . The results of the NNKI model within 
the SF region pertaining to the hysteretic properties, dy- 
namic symmetry breaking, and dynamic phase transi- 
tions (DPT) are summarized in the review article Ref. [5l 

In this paper we extend the NNKI model to include the 
next-nearest neighbor (NNN) interaction. We term this 
model the next-nearest neighbor kinetic Ising (NNNKI) 
model (see Eq. [T]). We study the effects of NNN in- 
teraction on the DPT (see Figs, [l]- [3|, correlation (see 



Fig. |4(a) ) , hysteresis loop area law and scaling (see 
Figs. [5] - 7|), metastable lifetime (see Figs, [s] - [9|, and 



order parameter fluctuation (see Fig. 10). There are two 
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main motivations behind this investigation. First, is the 
hysteresis loop area law effected by NNN interactions? 
Second, is the DPT boundary of the NNKI model sensi- 
tive to the NNN interactions? 

To obtain an expression for the hysteresis scaling law 
and to observe the DPT in the NNNKI model we study 
this model in the low frequency range (f G [10~^, 10~^]) 
at a constant temperature T=0.2Tc, where Tc is the crit- 
ical temperature of the zero field equilibrium NNN Ising 
model [24]. The physical justification for the choice of 
the frequency range is provided in Section |IVB| From 
our study we find that including the NNN interaction 
changes the hysteresis loop area law from h^-^^f^-^^ (NN) 
to hO-i4±o.oif0.i3±o.oi (NNN) and shifts the DPT bound- 
ary (see Figs. 1(a) and [2|. The order of the phase tran- 
sition is second order and this is confirmed by tracking 
the probability density of Q. We also study the effects of 
lattice size (L) on the hysteresis loop area law. Further- 
more, we explore the scaled variance of the dynamic 
order parameter to extract information on the critical 
properties of the NNNKI model. We also compute the 
density plot for the correlation function, B, and the hys- 
teresis loop area. A, at a fixed frequency for a range of 
ho and T values. 

The reason behind the change in the loop area law and 
the DPT boundary can be physically understood by ex- 
ploring the effects of NNN on the metastable lifetime, r 
(see Fig. |6|. The metastable lifetime is defined as the 
time taken by the system to decay to a net zero mag- 
netized state (m(t) = 0) from a completely magnetized 
state (m(t) = 1) upon an instantaneous field reversal. 
The sensitive dependence of r on L, h^ and T has been 
investigated thoroughly [25 . In this paper we focus 
on the effect of additional interactions on r. One of our 
main findings is that the metastable lifetime is dependent 
not only on ho, T, and L but also on additional interac- 
tions present in the model. This change is refiected in 
the value of the hysteretic exponents and in the location 
of the DPT phase boundary. Including the interactions 
however does not change the second order nature of the 
phase transition. Also, droplet theory as a mechanism 
for metastable decay still holds in the MD and SD region 
of the NNNKI model only. 

This paper is arranged as follows. In Se c. [TT| we pro- 
vide a review of the NNKI model. In Sec. IIIII we state 
the NNNKI model and describe the Monte Carlo method 
used to study its properties. In Sec. |IV| we present and 
discuss our results on NNNKI model. Finally, in Sec. [V| 
we summarize and state the main conclusions of our pa- 
per. 



II. NNKI REVIEW 

Unlike thermally quenched Ising systems the NNKI 
model is explcitily time-dependent [26l [27] . Hysteresis 



in this model arises out of a competition, R= 



_ 27r/a; 



be- 



tween two time scales [6] . The quantity R can be thought 
of as a scaled time or 1/R as a scaled frequency. Here 
2i\ jbj is the external field's period and r is the system's 
metastable life. The model's oscillatory nature causes the 
system to continuously alternate between momentarily 
stable and unstable states. The metastable state arises 
when the spins are aligned opposite to the field. When 
R< 1 the magnetization does not switch signs within a 
single period and the system has a non-zero value of Q 
with asymmetric hysteresis loops. However, if R>1 a 
symmetric hysteresis loop is observed and Q=0. 

Metastability has been utilized to show that there ex- 
ists two decay-mode regimes in the NNKI model - deter- 
ministic and stochastic ^Ij. The deterministic regime 
is split into the strong field (SF) region and the mul- 
tidroplet (MD) region. The stochastic regime contains 
the single droplet (SD) region and the coexistence (CE) 
region. Within either regime h^, L, and T not only de- 
termines the value of r but also plays a significant role in 
determining how the system decays from its metastable 
state. 

Droplet theory provides a physical explanation for the 
MD and SD regions. Systems within the MD or SD re- 
gions experience nucleation and growth of droplets of sta- 
ble spins over the entire metastable lifetime. Larger sys- 
tems experiencing MD decay will have many droplets cre- 
ated as it achieves stability. "Avrami's Law" provides a 
classical description of the MD region's mechanics [6][T2] . 
The SF region is characterized by "droplets" which are 
as large as the system size. In this region droplet the- 
ory breaks down. The stochastic regime is character- 
ized by weak external fields, small system sizes, and ex- 
tremely long average metastable lifetimes {{r{ho)) > 10^) 
[28 . The system within this region experiences a slow 
metastable decay as a single droplet of spin nucleates 
and eventually grows overtaking the system. 

The MD and SD regions have been studied exhaus- 
tively for finite size scaling effects, DPT, universality, and 
critical exponents [121 ESI IE] • It has been found that 
the critical exponents of the NNKI model are close to the 
equilibrium 2D Ising model and the 2D random percola- 
tion model [13]. The dynamical response of the NNKI 
model has also been studied subject to a square-wave 
oscillating external field with a soft Glauber dynamics 

[laEO]. 

In Ref. [5l the NNKI models hysteretic response has 
been studied using dynamic symmetry breaking. The 
hysteresis loop becoming asymmetric around the origin 
was used as a signature of DPT [5 . Rikvold and collab- 
orators used the metastability mechanism to show that a 
DPT within the NNKI model exists only when crossing 
between the MD and SF regions. Recently, evidence for 
a dynamic phase transition has been investigated exper- 
imentally in [Co/Ptja magnetic multilayers [31] . 
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III. MODEL AND METHOD 

The model used in this study is the kinetic NNN Ising 
ferromagnet on a square lattice with periodic boundary 
conditions. The NNNKI Hamiltonian is given by 

H = —Jnn ^ SiSj — Jnnn ^ SiSj — Hq sin(cjt) ^ Si 



(1) 



where is the ith spin and can have values of = ± 1, 
Jnn is the NN coupling, Jnnn is the NNN coupling, and 
ho is the external field amplitude. The sums Xl(ij) 

run over all the lattice sites, L^, for the NN and 
NNN interaction. Both the couplings are ferromagnetic 
and Jnn=Jnnn=l- The ratio of the couplings is defined 

as P=Jnnn/Jnn~l- 

The spin-flip dynamics used is the Metropolis algo- 
rithm with the Monte Carlo step per spin (MCSS) as the 
unit time step [32l |33l [SH [35]. The system is allowed 
to be in contact with a heat bath at temperature T, and 
each attempted spin flip from Si -Si is accepted with 
the probability W(Si -Si)=exp(-/3A Ei). Here /^Ei 
is the change in energy of the system that would result 
if the spin flip were accepted, and /3=l/kBT where the 
Boltzmann constant is set equal to one. 

Using the above Hamiltonian (Eq. [T]) and the Monte 
Carlo method we compute Q, the correlation (B), and 
the area (A) using 



2^ 



m{t)dt, 



ix) 



B = ^ (h m{t)h{t)dt, 

ZTT , 



A 



m{h)dh. 



(2) 
(3) 
(4) 



where m(t) is the time dependent magnetization per unit 
site. In the next two sections we study the physics of the 
NNKI model with NNN interactions. 



IV. RESULTS AND DISCUSSION 

In this section we report results from the Monte Carlo 
simulation of the NNNKI model. We calculate the DPT 
plot, P(Q) graphs, hystersis loop area exponents and scal- 
ing, correlation, and fluctuation of Q. The MCSS value 
used for these computations is 3.5x10^. Furthermore, in 
this paper we focus on the deterministic (SF/MD) region 
in the low frequency regime (f G 10~^ to 10~^). 



A. DPT and Probability distribution, P(Q) 



Fig. 1(a) shows the DPT plot for the NNNKI model. 
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FIG. 1: (Color online) Plots of the dynamic phase transition 
(DPT) and the probability density, P(Q), of the dynamic or- 
der parameter Q for L=64 at f=10~^. ho is the external field 
amplitude, T is the temperature, and the NN interaction is 
3nn=l- The letters A - I track the DPT along T=0.2Tc. 
(a) The overall qualitative feature of the DPT phase diagram 
for the NNN case is similar to the NN case (see Fig. 16c 
of Ref. 13 and Fig. |2] in this paper). The various phases: 
strong field (SF), multi-droplet (MD), and single droplet (SD) 
still survive. Quantitatively inclusion of the NNN interaction 
shifts the DPT phase boundary to higher values of magnetic 
field and temperature. The presence of additional interac- 
tions shift the DPT phase boundary and can give rise to an 
interaction induced DPT. (b) P(Q) plots for values of ho cor- 
responding to the letters A - I marked on the DPT plot. The 
X-axis represents Q and the y-axis P(Q). In the SD region 
there is a peak at Q—1 and in the SF region there is a peak 
at Q=0. In the intermediate MD region there is a double 
peak. The observation that the P(Q) displays no peak near 
Q=0 in the ordered dynamic phase is evidence that the DPT 
is a second order phase transition and it survives upto the 
NNN interaction. 



and T for a fixed frequency, f=10 ^. The figure dis- 
plays three regions - strong field ((|(5|)=0), multi-droplet 
(0< (IQI) <1), and single-droplet ((|Q|)=1). The cor- 
responding DPT plot for the NNKI model is shown in 
Fig. [2] Comparing the two plots we observe that while 
the overall shape is the same, the boundary of the DPT 
changes with the inclusion of NNN interactions. For the 
same choice of temperature the SD/MD region in the 
NNKI plot occupies a much smaller region compared to 



FIG. 2: (Color online) Dynamic phase transition plot for the 
NNKI model for L=64 and f=10~^. ho is the external field 
amplitude, T is the temperature, and the NN interaction is 
J^^ = l. Compared to the NNNKI DPT plot (see Fig. [l]) the 
SF/MD region occupies a much smaller space in the phase 
diagram. 
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FIG. 3: (Color online) Mean of the norm of the dynamic order 
parameter (|Q|) versus ho/Jnn corresponding to the letters A - 
I on the DPT plot of Fig. [l] ho is the external field amplitude 
and the NN interaction is Jnn=l- The values demonstrate 
the DPT from the SD to the MD and then finally to the SF 
region. The lines connecting the data points are a guide to 
the eye. 



the NNNKI-DPT plot. The system can undergo an in- 
teraction induced DPT. The sensitivity of the DPT phase 
boundary can be physically explained by considering the 
effect of interactions on the metastable lifetime of the 



system (see Section IV C). 

The nature of the DPT in the NNKI model was con- 
firmed to be second order by tracking the probablity den- 
sity of Q [12] . In this paper we follow the same procedure 
and compute the P(Q) graphs. The results are shown in 
Fig. [T(b)| There are nine plots which the track the evo- 
lution of P(Q) from the SD to the MD and then finally 
on to the SF region. Each plot is labelled by the letters 
A - I. Th e letters are also marked on the DPT plot (see 
Fig. |l(a)[ ). We see that deep in the SD region the prob- 
ability distribution has its weight shifted to Q = 1. As 
the value of ho is changed some of this weight is trans- 
ferred to values of Q which are not equal to one. This 
causes a bimodal distribution to be developed (plot (D) 




FIG. 4: (Color online) Density plots for the correlation, B, 
and the hysteresis loop area. A, at T=0.2Tc for L=64, f=10~^. 
The external field amplitude is given by ho, temperature by 
T, and the NN interaction by Jnn, where Jnn=l. (a) Density 
plot of the correlation. The system is maximally correlated 
in the SF region and minimally in the SD region, (b) Density 
plot for the hysteresis loop area. The system has the largest 
area for low temperatures and high magnetic fields. 



of Fig. l(b)[ ). This bimodal distribution rearranges itself 
furt her as we enter the MD region (see plots (E)-(F) of 
Fig. |l(b)'] ). Finally with a further increase in the value of 
ho the bimodal distribution weight shifts entirely to the 
Q=0 value with the peak centered around this point (see 
plots (G)-(I) of Fig.|l(b)D. This is in the SF region. 



The presence of the bimodal distribution of Q and the 
observation that P(Q) displays no peak near Q = in 
the ordered dynamic phase is a tell-tale sign of the sec- 
ond order phase transition present in the NNNKI system. 
We conclude from our computations that even though the 
additional interactions affect the phase boundary, the or- 
der of the phase transition is not affected by the presence 
of these interactions. 



In Fig. 4(a) we show the density plot for the correla- 
tion of the magnetization with the external field for the 
NNNKI model. From the plot we see that at large mag- 
netic fields and high temperatures the system is max- 
imally correlated. From the DPT plot the parameter 
regime in which this happens is the SF. The systems 
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(a)L=64 (a)Hysteresis plot, L=64,ho=16 
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(b)L=128 (b)Hysteresis plot, L=128,ho=16 



FIG. 5: (Color online) Log-Log plot of the frequency dis- 
persion of the hysteresis loop area law, A(ho,f), for the 
low frequency and the high magnetic field amplitude regime 
(ho > 8) of the NNNKI model at a temperature T=0.2Tc. f 
G [10~^,10~^] represents the frequency range and ho is the 
field amplitude. The choice of ho is such that the model is 
in the SF region. Inclusion of the NNN interaction causes a 
change in the hysteresis exponents. A(ho,f) is modified from 
(jv^jv^) j^o.i4±o.oif0.i3±o.oi (j^j^js^) ^he hysteresis 

loops corresponding to the area variation are shown in Fig. [4] 
The plots are shown for two lattice sizes. The lines connecting 
the data points are a guide to the eye. (a)L=64 and (b)L=128 



magnetization is least correlated with the external field 
for low fields and low temperatures. From the DPT graph 
this corresponds to the SD region. In Fig. 4(b) we show 
the density plot for the hysteresis loop area. The area has 
its maximum values in the SF region and its minimum 
values in the MD/SD region. At low temperatures and 
high field amplitudes the interaction energy dominates 
the system and the area is a maximum. However, as the 
temperature is increased thermal effects cause the loop 
area to diminish but remain finite. The SD region char- 
acterized by asymmetric hysteresis loops are the smallest 
in size and hence the area is small. 



B. Hysteresis 

In Fig.|5]we show the log-log plot of the frequency dis- 
persion of the hysteresis loop area law, A(ho,f), for the 



FIG. 6: (Color online) Hysteresis loops at T=0.2Tc and for 
ho=16 for various frequencies, f G [10~^, 10~^] represents the 
frequency range and ho is the field amplitude. The smallest 
loop area corresponds to f=10~^. The plots are shown for two 
lattice sizes. The dashed lines connecting the data points are 
a guide to the eye. (a)L=64 and (b)L=128 



low frequency range (f G [10 ^, 10 ^]) and the SF regime 
(ho > 8) of the NNNKI model. The results are shown 
for L=64 and L=128 at T=0.2Tc. For the NNKI model 
the low frequency range belongs to f <10~^. Above this 
range the NNKI model is in the high frequency regime 
with a different hysteresis loop area law [TB^ . Our Monte 
Carlo simulations of the NNNKI model over the fre- 
quency range f G [10-^, 10-^] has shown that for f> 10"^ 
the scaling law of the loop area breaks down. Further- 
more, hysteresis in the kinetic Ising model is dynamical 
in origin and disappears in the quasi-static limit. We 
therefore study the NNNKI model for f <10"^ [36^ and 
choose the frequency range f G [10~^, 10~^]. 

Using the Monte Carlo data we obtain the exponents 
in the low frequency range of the model. The results 
are summarized in Table [B We also show the standard 
deviation for each computed exponent. Since both the 
couplings are ferromagnetic, the system is aided by the 
presence of NNN interactions. The presence of NNN in- 
teractions in the model introduces additional competition 
for the external magnetic field to bring about a change. 
Therefore the response of the system to forming a hys- 
teresis loop of the same area as compared to the NN 
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FIG. 7: (Color online) Scaling laws for the loop area A(ho,f) 
at a constant temperature T=0.2Tc. The hysteresis area law 
data of Fig. [5] is replotted as function of the scaling variable 
j^o.i4f0.i3^ For the frequency range (f G [10-^ lO"^]) and the 
external field amplitude range ho =8 - 81 considered in this 
model the collapse of the data is good and the proposed area 
law is well obeyed, (a) Plot markers correspond to the var- 
ious values of ho=8, 16,32,48,64, and 81 (b)The plot markers 
correspond to the frequency range f G [lo-^lo-^]. 



TABLE I: Hysteresis exponents, a and b, for the loop area 
law, A(ho,f)(x hgf^ for the NNNKI model at a constant tem- 
perature. The exponents are computed from the Monte Carlo 
simulation results and then by a least-squares fitting of the 
data. The standard deviation for each computed exponent is 
shown. The results for the NN model are from Ref. 5 



Exponent 


NN model 


NNN model 


a 


0.70 


0.14 ± 0.01 


b 


0.36 


0.13 ± 0.01 



From Fig. [5] we observe that the loop area increases as 
either the frequency or the magnetic field increases. Such 
a behavior was observed for the NNKI model. Qualita- 
tively both the NNKI and the NNNKI loop areas display 
similar trends. In Fig. |6] we show the hysteresis loops 
corresponding to the external frequencies of Fig. |5] for 
L=64, 128 at ho=16. We see that as the frequency de- 
creases the loop area also decreases. Both Figs. 6(a)| and 
|6(b)| show similar shapes. In Fig. [t] we show the scaling 
of the hysteresis loop area. The two plots display the 
collapse of the data for either the field amplitude vari- 
atio n (see Fig. |7(a)[ ) or for the frequency variation (see 
Fig. 7(b)[ We observe that for the frequency range and 
field amplitudes considered the scaling relation for the 
loop area law is well obeyed. 



C. Metastable lifetime, t(/io, 7", L,p) 

The key to understanding the observed phenomena in 
the NNNKI model is to study the metastable lifetime. 
To determine r for this model we performed several in- 
stantaneous field reveral simulations [25]. The system 
was initially prepared in an all up spin configuration. 
The field was instantaneously reversed and the relax- 
ation of the system then studied. The metastable life- 
time is the number of MCSS needed for the system to 
decay to a net zero magnetized state {m{t) = 0) from 
a completely magnetized state {m{t) = ±1). The aver- 
age metastable lifetime, {r{ho^T^ L^p))^ is calculated af- 
ter 1000 repeated trials. The simulations were performed 
for L=64 at T=0.8Tc and ho=0.5Jnn for various ratios 
of the interaction strength p=0, 0.5,^,1. The results are 




model is diminished. This is refiected in the power of ho 
changing from 0.70 (NN) to 0.14 (NNN). The frequency 
exponent is similarly affected. Furthermore, to confirm 
that the computed exponents do not suffer from finite- 
size effects we have studied the lattice size dependence 
of the hysteresis exponents. We have computed the ex- 
ponents for L=50, 64, 72, 96, and 128. For these lattice 
sizes the a-exponents are all = 0.14 and the b-exponents 
are = 0.13. The values of the loop area law exponents in 
this paper are an average of these exponents. 



FIG. 8: (Color online) Variation of the average metastable 
lifetime, (r(/io, T, L,p)) as a function of interaction strength 
ratio, p=Jnnn/Jnn. The calculations were performed at 
T=0.8Tc and ho=0.5Jnn, where Jnn=l.We see from the plot 
that the value of r is sensitive to additional interactions 
present in the system. The physics of the NNKI and the 
NNNKI model is governed by the metastable lifetime. A 
change in the value of r causes the system to produce a shift 
in the DPT phase boundary and to obey a new hysteresis loop 
area law. The lines connecting the data points are a guide to 
the eye. 
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FIG. 9: Tracking the droplets for the NNNKI model at 
T=0.2Tc for L=200 using Metropolis dynamics. The sys- 
tem is initialized with all spins up, s^ = 1. The field was 
then instantaneously reversed and the relaxation of the sys- 
tem studied. The metastable lifetime is the number of MCSS 
needed for the system to decay to a net zero magnetized state 
(m(t) = 0) from a completely magnetized state (m(t) = ±1). 
Snapshots of the spin configurations in the SF, MD, and SD 
region are shown in the three rows respectively. The applied 
field was ho=-6Jnn (SF), -43 nn (MD), and -3.25 Jnn (SD). 
The first row represents the SF region at t=l, 2, and 3 MCSS. 
The second row represents MD at t=10, 20, and 30 MCSS. 
The third row represents SD at t=20, 30, and 40 MCSS. Sta- 
ble Si=-1 (metastable Si=-|-1) spins are represented by black 
(white). Droplet theory for the NNNKI model holds in the 
MD and the SD region only. 




FIG. 10: (Color online) Variation of the order parameter 
fluctuation x2=L^[{Q^) - (|Q|^)] versus ho/Jnn at f=10"^ 
and T=0.2Tc for L=64 where Jnn=l. The fluctuation 
pea ks aro und ho/Jnn=3.65 signifying the DPT transition (see 
Fig. 1(b)). The lines connecting the data points are a guide 
to the eye. 



the kinetic Ising model it is not obvious what the field 
conjugate to Q is. Recently, it has been demonstrated 
that the scaled variance of Q, X^, can be utilized as a 
proxy for the susceptibility in the kinetic Ising model to 
investigate the scaling behavior of the NEQ system near 
its critical period [37 . The variance of X is defined as 



L'm-{\Q\')] 



(5) 



In Fig. 



10 



we show the results for the fluctuation of Q, X^ 
with hoTJnn at T=0.2Tc and f=10-3 for L=64. The fluc- 
tuation peaks around ho/Jnn=3.65 signifying the DPT 



transition observed in Fig. |l(a) 



V. CONCLUSIONS 



shown in Figjsj 

From Fig. |8pwe see that the NNNKI model's average 
(r) is much greater than the NNKI model's metastable 
lifetime. Furthermore the lifetime is sensitive to the ra- 
tio of the interaction strengths. As mentioned earlier, 
hysteresis happens due to the competition between the 
metastable life of the system and the fields period. A 
change in the value of r will therefore cause the sys- 
tem to display a different DPT and hysteretic behavior. 
In Fig. |9] we show the formation of the droplets in the 
SD and MD region to confirm the decay mechanism of 
the NNNKI model. The SF region is not described by 
droplets. 

D. Order parameter fluctuation 

For equilibrium systems the fluctuation-dissipation 
theorem relates the susceptibility to fluctuations in the 
order parameter. As noted in earlier investigations for 



In this paper we explore the effects of NNN interac- 
tion in the kinetic Ising model. We study the effects on 
DPT and hysteresis loop area in the low frequency limit 
of the NNNKI model. We find that including the NNN 
interaction causes the DPT phase boundary to shift to 
greater values of magnetic field and temperatures. This 
can cause the system to undergo an interaction induced 
DPT. We also confirmed the second order nature of the 
DPT in the NNNKI model by studying the probabil- 
ity density P(Q) of the dynamic order paramter Q. We 
find that including the NNN interaction causes the hys- 
teresis loop area law to change from h^-^^f^-^^ (NN) to 
j^o.i4±o.oif0.i3±o.oi (NNN). One of our main findings is 
that the metastable lifetime, r, which dictates the physics 
of the NNNKI model is sensitive to additional interac- 
tions. Since the presence of additional interactions causes 
the value of r to change it causes the DPT phase bound- 
ary and the hysteretic properties to change. Finally, it 
is our hope that this work will motivate experimentalists 
to investigate the interaction induced DPT observed in 
the NNNKI model in real material systems. 
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